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Abstract 

We study the existence problem for a local implicit function determined by a system of nonlinear 
algebraic equations in the particular case when the determinant of its Jacobian matrix vanishes at 
the point under consideration. We present a system of sufficient conditions that implies existence of 
a local implicit function as well as another system of sufficient conditions that guarantees absence 
of a local implicit function. The results obtained are applied to proving new and classical results on 
flexibility and rigidity of polyhedra and frameworks. 
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1 Introduction 

Let F :Tl l x R m -> R" be a differentiable mapping, let t, i £ R 1 , let X, X E R m , and let F(t Q , X ) = 0. 
The classical Implicit Function Theorem provides conditions which imply that the equation F(t, X) = 
determines an implicit function X = X(t) in a neighborhood of the point (to, Xq). The principal condition 
is invertibility of the operator F' x (to, Xq). 

The Inverse Function Theorem has numerous applications and is generalized in various directions. 
However, the author is not aware of any version of this theorem which guarantees existence of an implicit 
function in the case when the operator F' x (to, Xq) is not invertible. In the present paper we will partially 
fill in this gap. 

Our study is motivated by that of flexible polyhedra and frameworks. It turns out that mappings 
F which appear is that field do not depend on the parameter t. We will focus our attention on this 
particular case. The following system of nonlinear algebraic equations can be considered as a typical 
example of a system to which our arguments can be applied: 

Fi(t,x\,x 2 ,x 3 ) = x\ + x\ - x\ - 1 = 0, 

F 2 (t, xi,x 2 ,x 3 ) = 3xi + x 2 - 3x 3 + 1 = 0, (1) 

F 3 (t,xi,x 2 , x 3 ) = xi - 3x 2 + x 3 + 3 = 0. 

The parameter t is not explicitly involved in this system. The point X — (5,5, 7) T satisfies (1). The 
determinant of the Jacobian matrix of (1) vanishes at X a = (x\, x, 2 , x 3 ): 



detF x (t,X ) 



2xi 2x 2 —2x 3 
3 1 -3 
1 -3 1 



10 10 -14 
3 1 -3 
1 -3 1 
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Therefore, the classical Implicit Function Theorem cannot be applied to (1). Nevertheless, from the 
results presented below it will follow that Xq is not an isolated solution to (1); on the contrary, it belongs 
to a continuous family of solutions X = X(t) which can be treated as an implicit function determined by 
system (1) and initial point Xq. 



2 Sufficient conditions for existence of an implicit function 

Let X = (xi, . . . , x m ) G R m and let F(X) = (-Fi(X), . . . , F n (X)), where each F k (k = 1, . . . ,n) is a 
polynomial. Without loss of generality, we may assume that the degree of each Fk is at most 2. 

To explain the last statement, we assume, for example, that in the system F(X) = under consider- 
ation each polynomial F k (k = 1, . . . , n — 1) is of degree at most 2 while F n (X) = x\xi — 1. Introduce 
a new independent variable x m +i and put X = {x\, . . . , x m ,x m +i). We also introduce new functions 
F n (X) = x m+1 x 2 -l and F n+1 = x m+1 - x\ and put F(X) = (F^X), . . . , F n _ x (X), F n (X), F n+1 (X)). 
Obviously, the system F(X) = is equivalent to the system F(X)) = and each equation of the latter 
system is of degree at most 2. 

Thus, without loss of generality, we may assume that the degree of each Ff. is at most 2. In this case, 
F k can written as 

mm m 

F k (X) = VV aZxiXj + V $Xi + 7 fe , 



i=l j=l 



i=l 



where , f3\ , and j k are some reals satisfying a*- 

It is well known that, if a system of polynomial equations admits a family of solutions which is con- 
tinuous with respect to a parameter, then this system also admits a family of solutions which depends 
analytically on a parameter (possibly different); see, for example, |l4| or Lemma 18.3 in (27) . Thus, assum- 
ing that the system F(X) = admits a continuous family of solutions X = X(t) = (xi(t), . . . , x m (t)) 7 we 
may assume without loss of generality that this family depends analytically on t, i.e., it can be expanded 
in a Maclaurin series: 



Xi{t) 



k=0 



x, k e R. 



Substituting this expansion into equation Fh(X) — 0, we obtain 



m m 

EE 

i=i j=i 



p=0 



f" 



OO 

q=0 



i=l 



ft 



p=0 



or 



p=0 



m m 

EE« 

i=l 3 = 1 



P 
q=0 



i" 



E EA fc 

p=0 L i=l 



+ 1 K 



Interpreting the left-hand side of the latter equation as the Maclaurin expansion of the function that 
equals zero identically, we conclude that, in this expansion, the coefficient of t p equals zero for each 
p > 1, i.e., the equation 



(2) 



E E E a % x i,i x 3,p-i + E ft* x *>p + 7 fc = o 

?'— 1 j=l <j— i— 1 

holds for all p > 1 and 1 < k < n. 

For each p > 1, put X p = X2, P , ■ ■ ■ , x mtP ) 6 R m . Let the bilinear mapping B : R m x R " — > R" 
be defined by the following rule: if X — (xi, . . . , x m ) £ R m and Y — (yi, . . . , y m ) G R m then the fcth 
coordinate of the vector B(X, Y) is equal to 



i=i j=i 
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Let the linear mapping A : R m — ► R™ be denned by the following rule: if X = (x\, . . . , x m ) G R m then 
the fcth coordinate of the vector A(X) equals 

m 
i=l 

Using this notation, we can rewrite (2) as 

£ B(X P , X q _ p ) + AX q = 0. 

p=0 

It follows that, if the vectors Xo, X\, . . . X q ^i are given and we seek the vector X g , then we need to 
solve the following system of linear equations: 

9-1 

B(Xq, X q ) + B(X q ,X ) + AX q = - B(X P , X q -p). (3) 

p=l 

Let the linear mapping C : R m ->■ R™ be given by CX = B(X ,X) + B(X, X ) + AX. Then we can 
rewrite (3) in concise form: 

9-1 

CX q = -J2B(X p ,X q _ p ). (4) 
P =i 

In the preceding consideration, the vectors X p were generated by the Maclaurin coefficients of a 
family of exact solutions to F(X) = 0. Now we assume that we have an arbitrary finite set Y , Y\, . . . , 
Y q of vectors in R m . We call the expression 

an approximate solution of degree q to the system of polynomial equations F{X) = if, for each p = 
1,2, ... ,q, the coefficient of t p in the Maclaurin expansion of the function F(Y(t)) is equal to zero. An 
equivalent formulation of this condition is as follows: for each p — l,2,...,q, the equation 

p-i 

CYp = -Y,B(Y l ,Yp_ l ) 
i=i 

holds. 

Now we are ready to formulate sufficient conditions implying existence of an implicit function which 
is determined by a system of algebraic polynomial equations. 
Theorem 1. Let 

(5) 

p=0 

be an approximate solution of degree q to a system of algebraic polynomial equations F(X) = 0. Suppose 
that there exists a number k (0 < k < q) such that, for all i = 1, 2, . . . , q and j = k, k + 1, . . . , q, the 
equation 

CY = -B(Y i ,Y j )-B(Y j ,Y i ) 
has a solution which lies in the linear span of the vectors Y^, Yfc+i, . . . Y q . Then F(X) = has an analytic 

oo 

family of solutions X(t) = ^2 X p t p whose initial coefficients coincide with the corresponding coefficients 
of the approximate solution (5), i.e., X p = Y p for each p = 0, 1, . . . , q. 
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Proof. Denote by L the linear span of the vectors Yfc, Yj 



Xn 



Let P be the set of all nonnegative integers p for each of which there exists an approximate solution 
of degree q + p, 

q+p 



(6) 



/=o 



to F(X) = such that (i) Xi = Yi for all I = 0, 1,. . . , q and (ii) Xi E L for each Z = g + 1, q+ 2, . . . , q + p. 

In view of the conditions of Theorem 1, G P. Hence, P ^ 0. Verify that P coincides with the set of 
all nonnegative integers, N. It suffices to show that if p E P then p + 1 G P. 

So, let p G P and let the approximate solution (6) satisfy (i) and (ii). To prove that p + 1 G P, it is 
sufficient to find a vector X q+p+ i G L satisfying the following system of linear algebraic equations: 



CX, 



q+p+l 



q+p 

£ 

1=1 



B{Xl,Xq +p+ \-l). 



(7) 



According to (i) and (ii), each of the vectors X q+1 , . . . ,X q+p lies in L and, hence, each of the vectors 
-Xg+i, ■ ■ ■ , X q+p is a linear combination of the vectors Yj., Yfc +1 , . . ., Y q . Therefore, the right-hand side of 
(7) is a linear combination of the vectors B(Yi, Yj) + B(Yj, Yi) with 1 < i < q and k < j < q. Then the 
conditions of Theorem 1 imply that there exists a solution to (7) which lies in L. Hence, p+leP and 
P = N. 

Thus, we see that F(X) = has approximate solutions of arbitrarily high degree whose initial coef- 
ficients coincide with the corresponding coefficients of the approximate solution (5). It remains to prove 
that these approximate solutions of arbitrarily high degree can be used to construct an exact solution in 
the form of a power series whose initial coefficients coincide with the corresponding coefficients of (5). 

Our approach is based on the following algebraic theorem by M. Artin (see || and |2(|): Given a 
system of polynomial equations f(x, y) = 0, where f = (/ x , . . . , f k ), x = (x\, . . . , x m ), y = (yi,..., y n ), 
there exists an integer (3 = /3(m,n, d, a) (depending on m, n, on the total degree d of the polynomials 
f, and on a nonnegative integer a) such that if f(x,y(x)) = (modx^), f3 = f3(m,n,d,a), for some 
polynomial y(x) then the system f{x,y) = has a solution y(x) that can be represented in the form of a 
convergent power series whose coefficients coincide with those of the polynomial y{x) up to the term x a . 

To complete the proof of Theorem 1, we apply Artin's theorem to F(X) — in the following way. 
Put a = q and let (3 be an integer whose existence is provided by Artin's theorem. From the above it 
follows that the approximate solution (5) can be extended to an approximate solution of an arbitrarily 
high degree, in particular, to an approximate solution of degree (3. Now Theorem 1 directly follows from 
Artin's theorem. 

We now discuss several examples of using Theorem 1. 

Example 1. Let F : R 3 -> R 3 be given by (1), namely, let 



Fi(t,x 1 ,x 2 ,x 3 ) =x\+xl 



1 = 0, 



F 2 (t,Xi, X2,X 3 ) = 3.Ti + x 2 - 3x 3 + 1 = 0, 

F 3 (t, xi,x 2 ,x 3 ) =xi — 3x 2 + X3 + 3 = Q 



and let Xq = (5, 5, 7) T . Direct calculations show that 



(4-) 







(4) = I 





(A 2 



7 2 = i; 
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(a?-)=^0 Oj; G9f)=^-3j; 7 3 

Y) = (x iyi + x 2 y 2 - x 3 y 3 , 0, 0) T ; 

\ / 10 10 -14 

1 -3 I; C= [ 3 1 -3 | : dot f = 

1 -3 1 

Solving the homogeneous system of linear algebraic equations CX = 0, we find that the vector X\ = 
(4, 3,5) T constitutes a basis for the space of its solutions. Direct calculations show that B(X\,X\) = 
(0,0, 0) T . Hence, we can put X q = for all q > 2. Thus, Theorem 1 can be applied with q = 2 and 
k = 1. It follows that Xq is not an isolated solution to F(X) = 0. In the case under consideration, the 
corresponding family of solutions may be explicitly written as X(t) = Xq + tX\. Its geometrical sense 
becomes obvious if we observe that the equation F\(X) = determines a one-sheet hyperboloid while the 
pair of linear equations F 2 (X) = F 3 (X) = determines its straight line generator which passes through 
the point Xq. 

Example 2. Let / : R 2 -> R 1 be given by the formula f(x\, x 2 ) — x\ — x\ and let Xq = (0, 0) T . 
Transform the equation f(X) = to a system of equations each of which is of degree 2: 

F 1 (x 1 ,x 2 ,x 3 ) = X\x 3 - x\ = 0, , , 

F 2 (xi,X2,X 3 ) = x\ - x 3 = 0. 



Then we have B(X,Y) = (^x x y 3 - x 2 y 2 + ±x 3 y-i, a;ij/i) 

C = 




0-1 

Instead of solving the corresponding equations (4) step by step, we note that the equation f(X) = has 
an obvious analytic family of solutions, x\ = t 2 , x 2 = t 3 . Whence we immediately obtain Xq = X\ = 
X 5 =X 6 = ... = (0, 0, 0)T, X 2 = (1, 0, 0) T , X 3 = (0, 1, 0) T , Xi = (0, 0, 1) T . Find the smallest values of 
q and k for which the hypotheses of Theorem 1 are fulfilled. 
Direct calculations show that 

B(X u Xi) +B(X i ,X 1 ) = (0,0) T for all i > 1, 

B(X 2 ,Xi) + B(Xi, X 2 ) = (0, 0) T for i = 3 or i > 5, 

B(X 3 ,Xi) + B{X U X 3 ) = (0,0) T for all % > 4, 

B(X 4 , Xi) + B(Xi, X 4 ) = (0, 0) T for all i > 4, 

B{X 2 ,X 2 ) = (0,1) T , B(X 3 ,X 3 ) = (-1,0) T , 

B{X 2 ,X 4 )+B{X^X 2 ) = (1,0) T . 

Hence, the hypotheses of Theorem 1 are fulfilled for q = k — 5 and are not fulfilled for any smaller values 
of q and k. Thus, as soon as we obtain the approximate solution X + tX\ + t 2 X 2 + t 3 X 3 + t 4 X 4 + t 5 X 5 , 
we can assert that it can be extended to an exact solution to (8). However, we cannot make the same 
assertion by using only the approximate solution X n + tX\ + t 2 X 2 + t 3 X 3 + t 4 X 4 . 

Now we give an example of a system of algebraic equations possessing an analytic family of solutions 
that cannot be obtained by Theorem 1 for any values of q and k. 

Example 3. Let F : R 3 R 2 b e given by the formulas 

F 1 (x 1 ,x 2 ,x 3 ) =x\ +x\+x\-A, 

F 2 { Xl ,x 2 ,x 3 ) = (xi - l) 2 + x\ - 1 U 
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and let Xq = (2, 0, 0) T . Direct calculations show that 

(#)= I I ; - l =-k 





^22/2 + X 3 y 3 , 35x2/1 + ^2y2) T 5 
~_ / 4 

° _ ^ 2 

Note that rankC = 1, imC = {(£,7?) G R 2 |£ = 2?7}, kerC = w) G R a |w = 0}, and dim kerC = 

2. 

Obviously, (9) determines the Viviani curve and thus admits the following analytic family of solutions: 

xi(t) = 1 + cosi, 
£2 (t) = skit, 
ar 3 (*) =2sm(£/2), 

OO 

X{t) = ( Xl (t),X2(t),x 3 (t)) = £ t p X p . 

N 

It is clear that each expression £ t p X p is an approximate solution of some degree to (9). Suppose it 

satisfies the hypotheses of Theorem 1 with some q and k. 

From the proof of Theorem 1 it follows that, for p > k, the vector X p is constructed as a linear 
combination of solutions to the equations 

CX = —B(Xi, Xj) — B(Xj,Xi) 

with 1 < i < q and k < j < q. Hence, each vector B(X2,X p ) + B(X p ,X2) must lie in the image of 
the operator C, i.e., its first coordinate must be twice as large as its second coordinate. However, this 
condition is not satisfied, since B(X2, Y) + B(Y, X2) — (— j/i, — j/i) T and infinitely many X p have nonzero 
first coordinate. 

Thus, example 3 shows that the conditions of Theorem 1 are not necessary for existence of an implicit 
function. This means that Theorem 1 cannot be used for proving that a given solution to a system of 
algebraic equations is isolated. In the next section we give several additional conditions under which the 
conditions of Theorem 1 are not only sufficient but also necessary for existence of an analytic family of 
solutions to a system of algebraic equations. 



3 Necessary conditions for existence of an implicit function 

A primary necessary condition for existence of a continuous family of solutions is known in the theory of 
bending of smooth surfaces at least since S. Cohn-Vossen ||. For systems of algebraic equations it can 
be formulated as follows (we use the notation introduced in the previous section): 

Theorem 2. If zero is the only solution to the system CX = then the system of algebraic equations 
F(X) = has no nonconstant analytic family of solutions representable in the form of a convergent power 
series with prescribed initial term Xq . 

Proof is carried out by way of contradiction. Suppose that F(X) = has a nonconstant analytic 
family of solutions which is represented in the form of a convergent power series 

00 

*(«) = JVXp 

p=0 
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and suppose that q is the smallest positive number such that X q / 0. According to (4), X q satisfies the 
following system of linear algebraic equations: 

9-1 

CX q = — ^ B(X p , X q ^p) = 0. 
p=l 

However, due to the hypotheses of Theorem 2, the last system has only zero solution. Hence, X q = 0. 
This contradiction proves Theorem 2. 

In the theory of bending of smooth surfaces some more advanced necessary conditions are also known 
(see, for example, B, An algebraic version is stated in the following theorem: 

Theorem 3. // the system of algebraic equations F{X) = and the vector Xq are such that no 
approximate solution of the first degree Xo+tXi, with X\ ^ 0, can be extended to an approximate solution 
of the second degree, then F(X) — has no nonconstant analytic family of solutions representable in the 
form of a convergent power series with initial term Xq- 

Proof is carried out by way of contradiction. Suppose that F(X) = has a nonconstant analytic 
family of solutions which is represented in the form of a convergent power series 

oo 

and suppose that q is the smallest positive number such that X q ^ 0. Then X q lies in the kernel of C 
and, according to the hypotheses of Theorem 3, it follows that B(X q , X q ) does not lie in the kernel of C . 
In view of (4), X 2q satisfies the following system of linear equations: 

2q-l 

CX2 q — — ^ B(X p , X 2q - P ) = —B(X q ,X q ). 
P =i 

Since B(X q ,X q ) does not lie in the kernel of C, the last system has no solutions. This contradiction 
proves Theorem 3. 

We now discuss several additional conditions under which the conditions of Theorem 1 are not only 
sufficient but also necessary for existence of an analytic family of solutions to a system of algebraic 
equations. The results we aim to obtain will generalize Theorems 2 and 3 in a sense. We start with the 
case in which there are few linearly independent vectors in the sequence X\, X 2 , ■ ■ ■ , X q of coefficients of 
approximate solutions. 

Theorem 4. Suppose that the system of algebraic equations F(X) — has an analytic family of 
solutions which is represented in the form of a convergent power series 

oo 

x(t) = J2t p x P 

and suppose that the vectors A3 and A4 belong to the linear span of the vectors X\ and X 2 ■ Then, for 
all 1 < i, j < 2, the system of linear equations 

CX — —B{Xi, Xj) — B(Xj,Xi) 

has a solution that belongs to the linear span of the vectors X\ and X 2 ■ 
Proof. The vectors Xi, X 2 , A3, and A4 satisfy the following equations: 

CAi = 0, 

CX 2 = -B(X 1 ,X 1 ), 

CX 3 = -B(X U X 2 ) - B(X 2 ,X 1 ), 

CXi = -B(X U X 3 ) - B{X 2 ,X 2 ) - B{X 3 , Ai). 
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Let L denote the linear span of X\ and X 2 . The second equation implies B(X\,X\) G CL. Since 
X 3 G L, the third equation implies B(Xi,X2) + B(X2,Xi) G CL. On the other hand, X 3 G L implies 
that X 3 = c\X\ + c 2 X 2 , with some reals c\ and c\. Hence, the fourth equation can be written as 
CX 4 = -24S(Xi,Xi) - cl[B{X u X 2 ) + B(X 2 ,X 1 )} - B(X 2 ,X 2 ). Here CX 4 belongs to CL according 
to the hypotheses of Theorem 4, B(X\,X\) and B{X ll X 2 ) + B(X 2 ,X\) belong to CL in view of the 
above proof. Therefore, B(X 2 ,X 2 ) G CL. This completes the proof of Theorem 4. 

Theorem 5. Suppose that the system of algebraic equations F(X) — has an analytic family of 
solutions which is represented in the form of a convergent power series 

00 

x(t) = J2t p x p 

p=0 

and suppose that the vectors X4, X 5 , X 6 , and X? belong to the linear span of the vectors X\, X 2 , and 
X 3 . Then, for all 1 < i,j < 3, the system of linear equations 

CX = -B(Xi,Xj) - B(Xj,Xi) 

has a solution that belongs to the linear span of the vectors X\, X 2 , and X3. 

Proof. Let a G R be an arbitrary real. Change the variable t = t + ar 2 in the solution X(t) = 

00 oc 

t p X p : Y(t) = X(t + oit 2 ) = T p Yp- Clearly, Y(t) is an analytic family of solutions to the equation 
F(Y) = and, thus, the equality 

9-1 

cy, = -^Bft,y„) (io) 
P =i 

holds for each q > 1. 

On the other hand, Y p can be written in terms of Xi by collecting similar terms in the expression 

00 00 

p=0 p=0 

We thus obtain 

Y = X Q ; 

Y 1= X i; (11) 

Y 2 = X 2 + aX u (12) 

Y 3 = X 3 + 2aX 2 ; (13) 

Y 4 = X 4 + 3aX 3 + a 2 X 2 ; (14) 

Y 5 =X 5 +4aX 4 + 3a 2 X 3 ; (15) 

Y 6 =X 6 + 5aX 5 + 6a 2 X 4 + a 3 X 3 . (16) 

According to the hypotheses of Theorem 5, each of the vectors X4, X5, and X§ belongs to the linear 
span of the vectors X\, X 2 , and X 3 . It follows that Xj = cjXi + c 2 X 2 + c 3 X 3 for each 4 < j < 6, with 
some reals c* (1 < i < 3, 4 < j < 6). Taking these formulas into account, we can write (14)~(16) in the 
following form: 

F 4 = (4 + 3")^3 + (cl + a 2 )X 2 + c\X 1} 

Y 5 = (c| + 4acl + 3a 2 )X 3 + (c 2 + 4acj)X 2 + (cl + 4acl)X u 

Y 6 = (c| + 5ac| + 6a 2 cl + a 3 )X 3 + (c 2 + 5ac 2 + 6a 2 4)X 2 [ '> 

+ (cl + bac\ + Qa 2 c\)X x . 



Let L denote the linear span of X\, X 2 , and X 3 . 
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For the case q = 2 equation (10) gives CY 2 = —B(Y\,Y\). Taking into account (11) and (12), we 
conclude from here that CX 2 + otCX x = -B(X 1 ,X 1 ), and finally that B(X 1 ,X 1 ) G CL. 

For the case q — 3 equation (10) gives CY 3 = —B(Y\,Y 2 ) — B{Y 2 ,Y\). Taking into account (12) and 
(13), we conclude from here that CX 3 + 2aCX 2 = -2aB(X 1 ,X 1 )-[B(X 1 ,X 2 ) + B(X 2 ,X 1 )], and finally 
that B(X 1 ,X 2 )+B(X 2l X 1 ) e CL. 

For the case q = 4, equation (10) yields CY A = -B(Y U Y 3 ) - B(Y 2 ,Y 2 ) - B(Y 3 ,Y 1 ). Taking (14) 
and (17) into account, we conclude that [c\ + 3a)CX 3 + {c\ + a 2 )CX 2 + c\CX x = -a 2 B(X 1 ,X 1 ) - 
3a[ J B(X 1 ,X 2 ) + J B(X 2 ,X 1 )] - [B{X U X 3 ) + B(X 3 ,X 1 )\ - B(X 2 ,X 2 ) and, finally, 

[B(X U X 2 ) + B(X 2 ,X 1 )] + B(X 2 ,X 2 ) G CL. (18) 

Similarly, from (10), we obtain in the case q — 5, 

[B(X 1 ,X 3 ) + B(X 3 ,X 1 )] + [B{X 2 ,X 3 ) + B(X 3 ,X 2 )} e CL, (19) 

in the case q = 6, 

(eg + 4acl + 3a 2 )[B{X u X 3 ) + B{X 3 , X x )\ 

+ {2c 2 - 10a 2 - Aacl)B(X 2 , X 2 ) (20) 
+ (cl + 3a)[B(X 2 ,X 3 ) + B(X 3 ,X 2 )} eCL, 

and, in the case q = 7, 

(ci - ac| + 6a 4 + c| + 5ac 3 5 + %a 2 4)[B{X u X 3 ) + B{X 3 ,X X )] 
+ (cl + 4acl - 2ac| + 8a 2 cl - 6a 3 )B(X 2 , X 2 ) 

+ (4 + 2acl + cl + 4a 2 )[B(X 2 ,X 3 ) + B(X 3 ,X 2 )} [ > 

+ {cl + 3a)B{X 3l X 3 ) e CL. 

We may treat relations (18)-(21) as a system of algebraic linear equations with respect to the following 
four vector-valued variables: B(X U X 3 )+B(X 3 , X^, B(X 2 ,X 2 ), B(X 2 , X 3 )+B(X 3 , X 2 ), &nAB{X 3n X 3 ). 
The right-hand sides of the corresponding equations are some vectors in CL. Denote them by U\, U 2 , 
U 3 , and U4. It suffices to prove that the determinant of the system is not equal to zero: in this case, 
each of the four vector- valued variables may be represented as a linear combination of U\, U 2l U 3 , and 
U4 and, thus, lies in CL. 

Direct calculations show that the determinant of the system corresponding to (18)— (21) equals 

-6a 4 + 18a 3 + 20a 2 c| + a[2c| + (c|) 2 - c|] 
+[-ci + c| - c 2 c| - {el) 2 + 2ctci - c|]. 

Obviously, this polynomial in a cannot be equal to zero identically for any values of the coefficients c*-. 
Hence, we can find a value of a such that the determinant of the system corresponding to (18)-(21) do 
not vanish. This completes the proof of Theorem 5. 

The above example 3 shows that, if the four vectors X\, X 2 , X 3 , and X4 are linearly independent, 
it may occur that some of the vectors B{X il Xj) + B{Xj 1 X i ) do not lie in the image of C while the 
system F(X) = determines some nonconstant implicit function. This means that there is no direct 
generalization of Theorems 4 and 5 to the case in which the four vectors X\, X 2l X 3 , and X4 are linearly 
independent. 

Obviously, this circumstance hampers proving that a system F{X) ~ admits no analytic family of 
solutions. There are also some other obstacles in proving this. The first obstacle is of purely technical 
nature and consists in a considerable increase of calculations: if we find that a first-order approximate 
solution Xq + tX\ admits an extension Xq + tX\ + t 2 X 2 to some second-order approximate solution 
then, for every X G kcr C, the expression Xq + tX\ + t 2 (X 2 + X) also gives us some second-order 
approximate solution; hence, we need to study the possibility of extending a second-order approximate 
solution to a third-order approximate solution for some family of second-order solutions rather than for 
a single solution. Another reason is of more principle character. Suppose that we are able to make our 
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way through the above-described increase of calculations and suppose we find a number TV such that no 
first-order approximate solution X + tX\, X\ e kcr C, X\ ^ 0, can be extended to any approximate 
solution of order N. Can we conclude that the system F(X) — defines no implicit function in a 
neighborhood of the point X n ? No, we cannot! We have to verify that there is no such an extension 
either for X + t ■ + t 2 X\, or for X + t ■ + t 2 ■ + t 3 Xi, or for any other approximate solution 
which has several vanishing initial coefficients (here Xi e ker C, X x ^ 0). Example 2 shows that a 
system F{X) = may admit an exact solution with several vanishing initial coefficients X p . On the 
other hand, we have no estimation for the number of zero coefficients in the Maclaurin expansion of an 
implicit function defined by a system F(X) = 0. Hence, we have to study an infinite set of cases caused 
by "writing zeros" at initial positions of an approximate solution. Therefore, in general, we have no 
algorithm which can guarantee absence of an implicit function. 

Nevertheless, below we will show that such an algorithm does exist in the case dim ker C = 1. First of 
all, we make the terminology more accurate. As before, let F(X) = be a system of algebraic equations 
each of which is of degree 1 or 2. Let a bilinear operator B and a linear operator C be constructed by 
means of the system. Suppose that dim ker C = 1. In the domain of C fix a codimension-1 subspace 

oo 

T such that T n ker C — {0}. A formal power series X(t) = t p X p is said to be a T -standard formal 

p=0 

solution to F(X) = if the following conditions arc fulfilled: 

9-1 

1) CX q = - £ B(X p ,X q _ p ) for each q > 1; 

P =i 

2) X\ ± 0; 

3) X p e T for each p > 2. 

The following theorem plays a key role in our approach: 

Theorem 6. Let the system F(X) = admit an exact nonconstant solution which can be represented 

oc 

in the form of a convergent power series X(t) = t p X P , let dim ker C — 1, and let T be a codimension-1 

oc 

subspace such that Tnker C = {0}. Then F(X) = admits a T -standard formal solution Y(t) = ^ t p Y p 

p=o 

such that Y = X . 

To avoid interrupting our presentation, we prove Theorem 6 at the end of this section. 

The coefficients Y p of a T-standard formal solution to F(X) = can be found as solutions to the 

p-i 

following system of linear algebraic equations: CY p = — Y^, B(Yi,Y p _i). The condition Y p G T implies 

i=i 

that the solution Y p is unique (if existent, of course) and thus no increase of calculations occurs. On 
the other hand, if, for some p, a solution Y p docs not exist then there is no T-standard formal solution 
to F(X) = 0. According to Theorem 6, this implies that F(X) = has no exact nonconstant solution 
in the form of a convergent power series (with arbitrarily many vanishing initial coefficients). Thus, we 
have a finite algorithm which, in some cases, can guarantee absence of an implicit function determined 
by F(X) = in a neighborhood of the point X - We present a test example of executing the algorithm 
proposed. 

Example 4. Let F : R 3 — ► R 3 be given by the formulas 

Fi(xi,x 2 ,x 3 ) = x\ + x\ + xl - 4, 
F 2 {xi,x 2 ,x 3 ) = (xi - 3) 2 +x 2 2 - 1, 
F 3 (xi,x 2 ,x 3 ) = x 2 

and let X = (2, 0, 0) T . It is clear that the equation F\ — defines a sphere in R 3 while F 2 = defines a 
cylinder which has a single common point with the sphere, namely, the point Xq. Thus F(X) — defines 
no implicit function. Demonstrate how Theorem 6 can be used to reach the same conclusion. 
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Direct calculations show that 
(««) 



(4) = 



(4) 

y) = (orij/i + a; 2 y2 + ^32/3, ariZ/i + x 2 y2, 0) T ; 

/ \ / 4 \ 

4= -6 J ; C = -2 . 
\ 1 / \ 1 / 

Note that rankC = 2, imC = {(£,»?, C) G R 3 |£ = 2ry}> ker C = {(«,", w) e R 3 |u = w = 0}, 
dim kcrC* = 1, X x = (0,0,1), and B(X 1 ,X 1 ) = (1,0, 0) T £ im C. 

The latter relation implies that the approximate solution X a + tX\ cannot be extended to any ap- 
proximate solution of the second order. Hence, there is no T-standard formal solution to F(X) = 0. By 
Theorem 6, F(X) = defines no implicit function in a neighborhood of Xq. 

Proof of Theorem 6. Let N be the least positive integer p such that X p ^ 0. Let q be the 

greatest integer such that the exact nonconstant solution to F(X) = (which, according to the claims 

00 

of Theorem 6, is representable in the form of a convergent power series X(t) = ^ t p X p ) possesses the 

p=0 

following properties: (i) X p = for all < p < q, p ^ (mod N) and (ii) X p 6 T for all < p < q, 
p =/= N, p = (mod N). 

We will verify that there exists a polynomial change of variables t = t(r) such that the new exact 

nonconstant solution X(t) = X(t(r)) to F(X) — 0, which is representable in the form of a convergent 

00 

power series X(t) = rP X pi possesses properties (i) and (ii) for q + 1 and is such that X = X , 

p=0 

Xn = Xn. 

In other words, we will prove that, by means of polynomial changes of variable t, the coefficients X p 
can be transformed one after another into the zero vector if p is not divisible by N and into some vectors 
lying in T if p is divisible by N in such a manner that X and Xn remain unchanged and, after such a 
transformation, we again obtain an exact solution to F(X) = representable in the form of a convergent 
power series. Accomplishing infinitely many such polynomial changes of variable t and not controlling 

the radii of convergence of power series which appear in this process, we obtain a formal power series 

00 

Z(t) = tP Z p whose coefficients Z p (p = 0, 1, . . .) possess the following properties: 

p=0 

(a) Z = X Q ; 

(b) Z p = for all p ^ (mod N); 

(c) Z N - X N + 0; 

(d) Z p GT for all p > 0, p = (mod N); 

9-1 

(e) CZ q = - B(Z P , Z q -p) for every q > 1. 

P =i 

Finally, executing the change of variable t = t n in the formal power series Z(t), we obtain a T- 
standard formal solution Y{t) whose existence is asserted in Theorem 6. 




(Pi) 



(Pi) 



(Pf) 




7 1 = -4; 



7 2 = -i; 



7 3 = 0; 
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So, to complete the proof of Theorem 6, it remains to prove that a solution X(t) possessing properties 
(i) and (ii) for some q, can be transformed into a solution X(t) possessing properties (i) and (ii) for q + 1 
and satisfying Xq = Xq, Xn = Xn- 

Put q = iN + j, where < j < N — 1. Consider a change of variable t = r + ar q+1 ~ N (here r is a 
new variable and a is a constant whose value will be specified later): 

X(t) =X(T + aT ( i+ 1 - N ) 

= X q + (t + ari+ 1 - N ) N X N + (r + OT<i+ 1 - N f N X 2N + -.. 

+ (r + a^+ 1 - N Y N X lN + (r + ar^")^ X iN+j +■■■ 
= X Q + t n X n + t 2N X 2N + ■■■ + T m X lN + T<HX q + NaX N ) + ■■■. 

First, consider the case < j < N — 1. We know that X p is a solution to the following system of 
linear algebraic equations: 

iN+j 

CX q = — ^ B(X r , X iN+ j_ r ). (22) 

r=l 

If r = (mod N) then iN + j — r = j (mod N) and, in particular, iN + j — r ^ (mod N). Hence, 
for arbitrary 1 < r < iN + j — 1, either X r = or Xjjv+j-r = 0. Therefore, the right-hand side of (22) 
equals zero and X q <G kerC. On the other hand, Xn G kcrC and dim kerC = 1. Consequently, the 
vectors X q and Xn are collinear. Since Xn ^ 0, there exists a such that X q + o:NXn = 0. Under such 
a choice of a, the exact solution X(t) = X(t + ar g+1 ~ N ) possesses properties (i) and (ii) for q + 1 as 
well as X = Xo and Xn = Xn- 

Now consider the case j = 0. In this case, if r = sN (0 < s < i) then iN + j — r = (i — s)N. Hence, 
(22) can be rewritten as 

i 

CX q = — B(X sN ,X^_ s ) N ). 

s=l 

Generally speaking, the right-hand side of the last expression does not equal zero. So, in general, X q 
does not lie in kcrC. Nevertheless, using the fact that the linear span of the subspaces T and kerC 
coincides with the whole space, we can find a such that X q + ciNXn G T. Under such a choice of a, the 
exact solution X(t) = X(t + ar q+1 ~ N ) possesses properties (i) and (ii) for q + 1 as well as Xo = Xo and 

Xn = Xn ■ _ 

Thus, the possibility has been proven of transferring the solution X(t) into a solution X(t) which 

possesses properties (i) and (ii) as well as X n = X 0} Xn = Xn- 
This completes the proof of Theorem 6. 

4 Applications to studying flexible polyhedra and frameworks 

Let if be a simplicial complex whose body is an (n — l)-dimensional connected compact topological 
manifold without boundary. A polyhedron in the n-dimensional Euclidean space R™ is, by definition, a 
continuous mapping / : K — > R™ which is linear on each simplex. Sometimes, the image of K under / is 
also referred to as a polyhedron. By a polyhedral sphere in R" we mean a polyhedron / : K — > R" with 
the body of K homomorphic to the sphere. 

We say that a polyhedron has no self-intersections if the mapping / is (globally) injective. In the 
present article, we consider polyhedra both with and without self-intersections. 

Definition. A polyhedron P = ,f{K) is flexible if there exists a family of polyhedra P t — (f t ,K), 
< t < 1, which is analytic with respect to the parameter t and for which the following conditions are 
satisfied: 

1) P = Po; 

2) for arbitrary 0-dimensional simplices Vj and Vk of K belonging to a 1-dimensional simplex of K, the 
equality \f(vj) — /(vfc)l = \ft( v j) ~ ft(vk)\ holds for all < t < 1 (henceforth \y\ stands for the Euclidean 
norm of a vector y = (yi, y 2 , y n ) e R™, i.e., \y\ 2 = y\ + y\ + . . . + y%); 
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3) there exist two O-dimensional simplices Vj and Vk of K which do not belong to any 1-dimensional 
simplex of K and for which the expression \ft(Vj) — ft(Vk)\ is n °t constant in t £ [0, 1]. 

A family P t with the above properties l)-3) is a nontrivial flexion of P. Note that the simplicial 
complex K remains constant during the flexion. 

In other words, a polyhedron is flexible if its spatial shape can be changed analytically with respect to 
a parameter (see condition 3)) without changing its intrinsic metrics (see condition 2)). Incidentally, the 
analyticity requirement with respect to the parameter can be substantially weekend. Namely, in flifl it is 
shown that if there exists a deformation of a polyhedron which is continuous with respect to a parameter 
and obeys the above conditions l)-3) then there also exists a deformation of the polyhedron which is 
analytic with respect to (possibly) another parameter and obeys conditions l)-3). 

During the last 25 years, the following two remarkable results were obtained in the theory of flexible 
polyhcdra: in 1977 R. Connelly gave an example of a flexible polyhedral sphere in R 3 without self- 
intersections 0, and in 1996 I. Kh. Sabitov published a complete proof of the statement asserting that 
each flexible polyhedron (even with self-intersections) in R 3 preserves the (oriented) volume bounded 
by it during the flexion p2]| . The last statement was known during several decades as the "Bellows 
conjecture". Other proofs of this statement may be found in |u|, |23]] , and p3 . 

We are interested in the problem of "practical" recognition whether a given polyhedron / : K — > R" 
is flexible or not. It is sufficient to describe a flexion /j : if — > R" of the polyhedron via motions of its 
O-dimensional simplices Vj (1 < j < N): xj(t) = ft(vj) £ R™. Furthermore, the above conditions l)-3) 
can be reformulated as follows: 1') the vectors Xj(0) are given; 2') if O-dimensional simplices Vi and Vj 
are joint together in K by a 1-dimensional simplex then the equality \xi(t) — Xj(t)\ 2 — \xi(0) — Xj(0)\ 2 
holds true for all < t < 1; 3') there exist O-dimensional simplices Vi and Vj that are not joint together 
in K by a 1-dimensional simplex and such that the expression \xi(t) — Xj(t)\ 2 is not constant in t £ [0, 1]. 

In other words, the problem of whether a given polyhedron is flexible or not is equivalent to the 
problem of whether the set of vectors xj = Xj(0), j = 1, 2, . . . , N, is an isolated solution to the system of 
algebraic equations 

\xi -x 3 \ 2 = \ Xi (0) -Xj(0)| 2 (23) 

or this system defines an implicit function xj = Xj(t) in a neighborhood of the point xj = Xj(0), 
7 = 1,2,..., N. It is worth bearing in mind that we are not interested in motions of f(K) as a rigid body 
in R", i.e., we are only looking for solutions with properties 3) or 3'). This requirement can be easily 
satisfied in the following manner. Fix an (n — l)-dimensional simplex of K . Suppose its O-dimensional 
simplices are denoted by v\, V2, ■ ■ . , v n . We agree that ft{v\) lies at the origin during the course of 
deformation of the polyhedron (and, thus, always has zero coordinates); ftfa) always lies on the first 
coordinate axis in R™ (and, thus, all but the first coordinates of it vanish identically); ft{v?,) always lies 
in the 2-dimensional plane spanned by the first and second coordinate axis in R" (and, thus, all but the 
first and second coordinates of it vanish identically) ; and so on. This construction reduces the number 
of independent variables in (23) but does not preclude the above reduction of the decision problem of 
whether the polyhedron is flexible to the problem of whether the given solution to the system of algebraic 
equations is isolated or the system defines an implicit function in a neighborhood of the solution. 

The decision problem of whether a framework is flexible in R™ can be reformulated in a similar way. 

A framework in R" is a connected graph whose vertices are points in R™ and edges are straight line 
segments joining some of its vertices. It is conventional to call the vertices of a framework joints and 
straight line segments bars. The set of 0- and 1-dimensional faces of a polyhedron in R™ can be considered 
as a typical example of a framework. 

A framework in R" is said to be flexible if it admits a nontrivial analytic deformation, i.e., the positions 
of its joints can be changed in R™ analytically with respect to a parameter in such a way that the length 
of each bar remains constant while the distance between some two joints (which are not joint together 
by a bar) is not constant. 

Sometimes, it is also useful to study so-called pinched frameworks, i. e. such frameworks that the 
spatial positions of some of their joints are fixed and should not be changed during deformations. 

As we have just mentioned above, a suitable framework can be associated with any polyhedron 
which consists of the set of its 0- and 1-dimensional faces. This framework is called a 1-skeleton of the 
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polyhedron. Obviously, a polyhedron is flexible if and only if its 1-skeleton is flexible (the reason is that, 
according to our definition, each face of a polyhedron is a simplex). Hence, the decision problem of 
whether a given polyhedron is flexible or not is a particular case of the decision problem of whether a 
given framework is flexible or not. So, we will focus our attention on the latter problem. 

The spatial position of a framework is determined as soon as the positions of its joints Xj(0) € R ra 
are given. The decision problem of whether a given framework is flexible or not is, obviously, equivalent 
to the problem of whether the set of vectors Xi(0) is an isolated solution to (23) or defines an implicit 
function in a neighborhood of 2^(0). Furthermore, is necessary to exclude trivial deformations when 
the framework moves as a rigid body. This can be done in the same way as it was previously done for 
polyhedra. This problem does not appear for pinched frameworks at all, since they often do not permit 
any trivial deformations. 

In Section 2, we introduced the notion of an approximate solution of degree q to a system of polynomial 
equations. As applied to system (23) associated with a framework, the term "infinitesimal flexion of order 
q" is conventionally used. More precisely, let a framework in R™ be determined by the positions Xj(0), 
i = 1,2, ... ,N, of its joints. The set of vectors Xi tP , i — 1, 2, . . . , N, p = 1, 2, . . . , q is said to be an 
infinitesimal flexion of order q of the framework if 



= (mod t q ) 



for all indices i and j such that the joints Xi(0) and Xj(0) are joint by a bar. An infinitesimal flexion 
is said to be trivial if it is an initial part of the Taylor expansion of the trajectories of Xi (0) under the 
action of some one-parameter group of isometries of R™. 

A framework is said to be qth-order infinitesimally flexible if it admits a nontrivial infinitesimal flexion 
of order q. Otherwise it is called qth-order infinitesimally rigid. 

Roughly speaking, the following theorem asserts that if a framework admits a "regular" infinitesimal 
flexion of sufficiently large order then it is flexible. 

Theorem 7. Let a framework P be qth-order infinitesimally flexible and let X]p=o X p t p be a nontrivial 
infinitesimal flexion of order q. Let operators B and C be constructed for system (23) which is associated 
with P and let there exist a number k (0 < k < q) such that, for every i = 1,2, ...,<? and every 
j = k, k + 1, . . . , q, the equation 

CX = —B(Xi, Xj) — B(Xj,Xi) 

has a solution lying in the linear span of the vectors Xk, Xk+i, ■ ■ ■ , X q . Then P is flexible. 
Proof of Theorem 7 follows immediately from Theorem 1. 

All flexible octahedra in R 3 were classified by R. Bricard [|| (see also [|o|). In 0, it is shown (in 
slightly different terms) that the conditions of Theorem 7 hold true with q — 5 and k = 1 for the 1-skeleton 
of so-called Bricard's flexible octahedra of the first type. 

The idea to use infinitesimal rigidity (of some order) of a framework for proving its rigidity is explored 
for a long time and is based, first of all, on the following 

Theorem 8. Every 1- order infinitesimally rigid framework in R™ is rigid. 

Proof follows immediately from Theorem 2 and from the above reduction of the decision problem of 
whether a framework is flexible to the problem of whether a given solution Xj — Xj(0), j — 1, . . . , N, to 
system (23) of algebraic equations is isolated or the system defines an implicit function in a neighborhood 
of the solution. 

Theorem 8 is one of the corner stones in problems of existence and uniqueness for convex polyhedra in 
the way of exposition which is used in the classical book B. Among recent papers that use Theorem 8, we 
indicate the articles [jl{| and [|l6| where it is shown that both in R 2 and R 3 there exist rigid triangle-free 
frameworks with all bars having length 1. Other statements about interrelations between rigidity and 
infinitesimal rigidity may be found in [|, Q, @, H§, §§, |o), H- 

The following theorem was proven for the first time (by different methods) in H|. A similar theorem 
for smooth surfaces was obtained in Q. 

Theorem 9. Every 2- order infinitesimally rigid framework in R™ is rigid. 
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Proof of Theorem 9 ensues immediately from Theorem 3. 
The following theorem generalizes Theorems 8 and 9. 

Theorem 10. Let a framework K in R™ have a single nontrivial linearly independent first-order 
infinitesimal flexion and let there exist a number q > 1 such that K is qth-order infinitesimally rigid. 
Then K is rigid. 

Proof. Eliminate trivial motions as was described above. The kernel of the operator C which is 
associated with (23) has dimension 1. Put T = (kerC) . According to Theorem 6, (23) has a T-standard 
formal solution Y(t) = J2p= Y p t p such that Y = K and K + Y x t is a nontrivial first-order infinitesimal 
flexion. Furthermore, Y(t) = ^2p =0 Y p t p is a nontrivial gth-order infinitesimal flexion. However, this 
contradicts the hypothesis of Theorem 10 according to which K is qth-order infinitesimally rigid. This 
contradiction proves Theorem 10. 

For smooth surfaces, a theorem similar to Theorem 10 was obtained (by different methods) in plj . 
For some other results about interrelations between higher-order infinitesimal rigidity and rigidity see, 
for example, p5fl , and Q (for frameworks) and Q and [|l9| (for smooth surfaces). 

References 

[1] Aleksandrov A D (1958) Konvexe Polyeder. Berlin: Akademie-Verlag 

[2] Alexandrov V (1998) Sufficient conditions for the extendibility of an n-th order flex of polyhedra. 
Beitr. Algebra Geom. 39: 367-378 

[3] Artin M (1969) Algebraic approximation of structures over complete local rings. Publ. Math., Inst. 
Hautes Etud. Sci. 36: 23-58 

[4] Bricard R (1897) Memoire sur la theorie de l'octaedre articule. J. Math. Pures Appl. 3: 113-148 

[5] Bolker E, Roth B (1980) When is a bipartite graph a rigid framework? Pacific J. Math. 90: 27-44 

[6] Cohn-Vossen S (1936) Bending of surfaces in the large. (Russian). Usp. Mat. Nauk 1: 33-76 

[7] Connelly R (1977) A counterexample to the rigidity conjecture for polyhedra. Publ. Math., Inst. 
Hautes Etud. Sci. 47: 333-338 

[8] Connelly R (1980) The rigidity of certain cabled frameworks and the second-order rigidity of arbi- 
trarily triangulated convex surfaces. Adv. Math. 37: 272-299 

[9] Connelly R, Whiteley W (1996) Second-order rigidity and prestress stability for tensegrity frame- 
works. SIAM J. Discrete Math. 9: 453-491 

[10] Connelly R, Sabitov I, Walz A (1997) The Bellows conjecture. Beitr. Algebra Geom. 38: 1-10 

[1 1] Connelly R, Servatius H (1994) Higher-order rigidity — what is the proper definition? Discr. Comput. 
Geom. 11: 193-200 

[12] Efimov N V (1948) Qualitative problems in the theory of bendings of surfaces. (Russian). Usp. Mat. 
Nauk 3, no.2: 47-158 

[13] Efimov N V (1952) Some statements about infinitesimal rigidity and rigidity. (Russian). Usp. Mat. 
Nauk 7, no.5: 215-224 

[14] Gluck H (1975) Almost all simply connected closed surfaces are rigid. Lect. Notes Math. 438: 225-239 

[15] Maehara H, Chincn K (1995) An infinitesimally rigid unit-bar-framcwork in the plane which contains 
no triangle. Ryukyu Math. J. 8: 37-41 



Implicit Function Theorem for Systems of Polynomial Equations 



18 



[16] Maehara H, Norihide T (1996) A spatial unit-bar-framework which is rigid and triangle-free. Graphs 
Comb. 12: 341-344 

[17] Maehara H (1998) Vector fields and quadratic surfaces. Ryukyu Math. J. 11: 53-63 

[18] Perlova N G (1991) On the interrelation between /cth-order rigidity and analytic inflexibility of 
surfaces. (Russian). Ukr. Gcom. Sb. 34: 98-104 

[19] Perlova N G (1995) On the interrelation between the rigidity of order k > 3 and analytic deformability 
of surfaces. (Russian). Mat. Fiz. Anal. Gcom. 2: 456-462 

[20] Sabitov I Kh (1992) Local theory on bondings of surfaces. Geometry III. Theory of surfaces. Encycl. 
Math. Sci. 48: 179-250 

[21] Sabitov I Kh (1994) On the relations between infinitesimal bendings of different orders. J. Math. 
Sci., New York 72: 3237-3241 

[22] Sabitov I Kh (1996) The volume of a polyhedron as a function of its metric. (Russian). Fundam. 
Prikl. Mat. 2: 1235-1246 

[23] Sabitov I Kh (1998) A generalized Heron- Tartaglia formula and some of its consequences. Sb. Math. 
189: 1533-1561 

[24] Sabitov I Kh (1998) The volume as a metric invariant of polyhedra. Discrete Comput. Geom. 20: 
405-425 

[25] Stachel H (1999) Higher-order flexibility for a bipartite planar framework. Preprint. Vienna Univ 
Technol 

[26] Stachel H (1999) Infinitesimal flexibility of higher order for a planar parallel manipulator. Preprint. 
Vienna Univ Technol 

[27] Wallace A H (1958) Algebraic approximation of curves. Can. J. Math. 10: 242-278 

[28] Whiteley W (1984) Infinitesimally rigid polyhedra. I: Statics of frameworks. Trans. Amcr. Math. 
Soc. 285: 431-465 

[29] Whiteley W (1984) Infinitesimal motions of a bipartite framework. Pacific J. Math. 110: 233-255 

[30] Whiteley W (1985) The projective geometry of rigid frameworks. In Batten L, Baker C (eds.) Finite 
Geometries. Marcel Dekker: New York, pp 353-370 

[31] Whiteley W (1987) Rigidity and polarity. I: Statics of sheet structures. Geom. Dcdicata 22: 329-362 



